THE ANISOTROPIC TRANSPORT EFFECTS ON THE 

DILUTE PLASMAS 



EBRU DEVLEN 

University of Ege, Faculty of Science, Department of Astronomy & Space Sciences, 

Bornova, 35100, Izmir, TURKEY 

ebru . devlenSege . edu . tr 
Received ; accepted 



- 2 - 
ABSTRACT 

We examine the linear stability analysis of a hot, dilute and differentially ro- 
tating plasma by considering anisotropic transport effects. In the dilute plasmas, 
the ion Larmor radius is small compared with its coUisional mean free path. In 
this case, the transport of heat and momentum along the magnetic field lines 
become important. This paper presents a novel linear instability that may more 
powerful and greater than ideal magnetothermal instability (MTI) and ideal mag- 
netorotational instability (MRI) in the dilute astrophysical plasmas. This type 
of plasma is believed to be found in the intracluster medium of galaxy clusters 
and radiatively ineffective accretion flows around black holes. We derive the dis- 
persion relation of this instability and obtain the instability condition. There 
is at least one unstable mode that is independent of the temperature gradient 
direction for a helical magnetic field geometry. This novel instability is driven 
by the gyroviscosity coupled with differential rotation. Therefore we call it as 
gyroviscous modified magnetorotational instability (GvMRI). We examine how 
the instability depends on signs of the temperature gradient and the gyroviscos- 
ity, and also on the magnitude of the thermal frequency and on the values of the 
pitch angle. We provide a detailed physical interpretation of obtained results. 
The GvMRI is apphcable not only to the accretion flows and intracluster medium 
but also to the transition region between cool dense gas and the hot low-density 
plasma in stellar coronae, accretion disks, and the multiphase interstellar medium 
because of being independent of the temperature gradient direction. 

Subject headings: accretion, accretion disks — MHD — dilute plasmas — gyroviscosity — 
thermal conduction — parallel viscosity-black holes 



-3- 



1. INTRODUCTION 

In recent years, MRI has been recognized as a powerful source of angular momentum 
transport in accretion discs. An accretion disc with angular velocity decreasing outward 
and threaded by a weak magnetic field is linearly unstable (Balbus & Hawley 1991). Local 
and global simulations of Keplerian discs showed that the MRI leads to a turbulence and 
therefore energy and angular momentum is transported outward (Balbus 2003). Similiarly, 
the dilute, stratified plasma is bouyantly unstable when the temperature increases in the 
direction of gravity. Heat is transported mainly along magnetic field lines in such a medium 
(Balbus 2000, 2001). The MTI has been studied with nonhnear simulations. Parrish 
and Stone (2005) investigated the nonhnear evolution of the MTI and showed that the 
instability causes turbulence and heat transport. They noted that MTI may explain the 
almost isothermal temperature profile observed in the outer part of X-ray emitting regions 
in intracluster medium (ICM) of galaxy clusters and the structure of radiatively inefficient 
accretion flows. In the MRI and MTI, the weak magnetic fields act as to turn free-energy 
gradients into sources of the instability. 

If the plasma is sufficiently dilute, the viscous stress tensor is also anisotropic 
(Braginskii 1965). Balbus (2004) showed that the viscous stress tensor can cause a strong 
instability in the dilute astrophysical discs. He showed that the maximum growth rate of 
instabihty exceeds that of the MRI. But, he only took parallel component of the stress 
tensor into account. The cause of the magnetoviscous instability is the same with MTI: 
Initially magnetic field fines are isorotational (isothermal in the case of MTI). Perturbed 
magnetic field lines are stretched out in the direction of the angular velocity (temperature) 
gradient. Thus, angular momentum (heat) is transferred from one fluid element to another. 
One fluid element at smaller radius drops down to smaller radii and other element moves to 
more distant radii. So, the field lines become more bent and the process runs away (Islam 
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& Balbus 2005). 

Indeed, to understand the true nature of the dilute astrophysical plasmas, we consider 
magnetohydrodynamic (MHD) equations which include the terms describing transport 
of heat and momentum by thermal conduction and viscosity (Braginskii 1965). Dilute 
means that the ion Larmor radius (r^i) is small compared with a mean free path (Aj) 
and any macroscopic length scales in the plasma. This is tantamount to saying that the 
ion cyclotron frequency (ujci) is much greater than the ion-ion collision frequency {ui ). 
Under this condition, parallel heat conduction of electrons is much larger than that of 
ions by the factor (mj/mg)^/^ and the parallel viscosity of ions is much larger than that 
of electrons by the same factor. Therefore, an anisotropic electron heat conduction and 
an anisotropic ion viscosity must be taken into account in the MHD equations. Ramos 
(2003) obtained dynamic evolution equations of the parallel heat fluxes in a coUisionless 
magnetized plasma. Besides, he noted that neglecting the parallel heat fluxes in low 
coUisional regime cannot be justified physically. He emphasized that one must consider 
the contribution of the gyroviscosity in the stress tensor for consistency in the analysis. 
Ramos (2005) presented the fluid moment equations with Finite Larmor Radius (FLR) 
effect for coUisionless magnetized plasmas. His analysis included the gyroviscous stress, 
the pressure anisotropy and the anisotropic heat fluxes. He claimed that his formalism 
can be applicable for arbitrary magnetic fleld geometry and arbitrary plasma pressure and 
also fully electromagnetic nonlinear dynamics. He extended previous study considering 
coUisional terms based on full Fokker-Planck operators for non-Maxwellian distribution 
functions. The low coUisional regime of interest is described with two small parameters: the 
ratio of the electron to the ion masses is comparable to S (i.e., {me/rrii)^^'^ ^ ^ ^ ^) ^-iid 
the ratio of the ion coUision to cyclotron frequencies is smaller than 5^ (i.e., Ui/ud < 5^), 
where 5 is the fundamental expansion parameter which is the ratio of the ion Larmor radius 
to the shortest macroscopic length scale (Ramos 2007). 
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Recently, Ferraro (2007) examined the FLR effects on tfie magnetorotational instability. 
He showed that the FLR effects are dominant compared to the other effects which includes 
the Hall effect in the limit of weak magnetic fields. He restricted his analysis in a vertical 
magnetic field geometry. He found that the growth rate of unstable mode is around 
when the ratio of the gyroviscous force to the magnetic tension force is greater than zero. 
But, when the ratio of the gyroviscous force to the magnetic tension force is smaller than 
zero, there is no unstable mode. Devlen & Pekiinlii (2010, Paper I) investigated the 
stability properties of weakly magnetized, dilute plasmas by considering combined effects of 
gyroviscosity and parallel viscosity which are components of the stress tensor in the presence 
of a helical magnetic field geometry. They showed that although the parallel viscosity is 
greater than gyroviscosity under the condition of dilute plasma, it hasn't any effect on 
the instability condition and growth rates. Also they showed that the powerful instability 
emerges due to finite Larmor Radius effects. They estimated that the growth rates of this 
GvMRI varied in the range of O.bfl — 3fl for the different values of pitch angles which is the 
angle between the magnetic field vector and the (f) axis of the coordinate system. When the 
ratio of the gyroviscous force to the magnetic tension force is smaller than zero (i.e., in the 
case of O B^), they found that there are unstable modes with growth rates around 2fl. 
This result was contrary to that of Ferraro's (2007). 

To clearly comprehend the dynamics of dilute astrophysical plasmas, all the anisotropic 
transport effects must be taken into account. Therefore, in this work, we extend our 
previous study (Paper I) by considering anisotropic electron heat conduction term. We show 
that weakly magnetized, differentially rotating dilute plasmas are unstable in the presence 
of the parallel viscosity, gyroviscosity and thermal conduction. This novel instability is 
extremely powerful and occurs at all wavenumbers. 

This paper is organized as follows. In the next section (§2), we give the linearization of 
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the MHD equations used in our analysis. We examine the physical structure of modes. We 
derive the dispersion relation of instability and the instability criterion. In §3, we examine 
the numerical solutions of the dimensionless dispersion relation and finally in §4, we discuss 
the physical interpretation of instability and summarize the our results . 



The plasma fulfilling the conditions r^j <^ Aj and e = UciTi ^ 1 is called dilute plasma, 
where is the ion Larmor radius, Aj is the ion collision mean free path, coci is the ion 
cyclotron frequency and Ti = l/vi is the inverse of the ion- ion collision frequency. The 
presence of the magnetic field introduces anisotropy to the medium. Cyclotron frequencies 
of the plasma species and the velocity gradients at macroscopic scales are the sources of 
anisotropy. 

If plasma consists only of hydrogens then e may be taken as (Spitzer 1962) 



where n is the proton number density in cm~^, T4 the temperature in units of 10^ K, 
Bfj^G is the magnetic field in units of 10'^ and In A is the Coulomb logarithm. The condition 
of e ^ 1 is fulfilled even in the presence of a very weak field with n <1 and T4 > 1. 

Under these conditions, the plasma dynamics described by MHD equations should 
include the anisotropic terms accounting for the free flow of particles along the magnetic 
field hues (Braginskii 1965). Ion parallel viscosity is higher by a factor (mj/mg)^''^ than 
that of electrons. So, the viscosity of the dilute plasma is determined mainly by the ions. 



2. 



LINEAR STABILITY ANALYSIS 



2.1. Dilute Plasma Properties 




(1) 
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Even if the ion viscosity is very small, in a rotating system, it may become very important 
(Balbus 2004). Similarly, since the electron contribution to the heat flux is higher than 
the ion contribution by a factor of (mj/me)^/^, the ion contribution may be considered as 
negligible. Since the electrons have mean free paths much longer than their gyro-radii in the 
dilute plasma, the thermal conductivity is strongly anisotropic. That is, in the astrophysical 
dilute plasma threaded even by a weak magnetic field, the momentum by ions and heat flux 
by electrons is transported primarily along the magnetic fleld lines. 



In an attempt to investigate parallel viscosity, gyroviscosity and heat flux in a dilute 
plasma, one should consider the two-fluid equations. Below are the standart extended MHD 
equations which arc obtained by using two-fluid equations including stress tensor, 11 , and 
the heat flux Q (see Appendix A): 



2.2. Basic Equations 




(2) 




(V X B) X B 



(3) 



c 



dB 

'dt 



V X (v X B) , 



(4) 



^ + ^P(V.v)^4v.Q, (5) 

where p is the mass density, v is the fluid velocity, P is the scalar pressure, 11 is the 
stress tensor, B is the magnetic fleld, g is the gravitational acceleration, Q is the heat flux. 
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and d/dt — d/dt + V ■ V is & Lagrangian derivative. 

Stress tensor have three components which are called as parallel (||), perpendicular 
(±) and the gyroviscous (gv) (Braginskii 1965). Perpendicular viscosity is greater than the 
parallel viscosity by a factor of {vLi/l)^, therefore it may not be taken into account; where 
ru is the Larmor radius and / is the mean free path of the particles in a dilute plasma. So, 
we used parallel and gyroviscous components of the stress tensor which are given by 

n" = 0.96^ (l - 3bb) (b • W • b) , (6) 



p. 

jjgv I 



b X W • (l + 3bb) + 



b X W • + 3bb) 



(7) 



where b = B/S, ujd = eB/rUiC are the unit vector along the magnetic field and the 
cyclotron frequency. Ui is the ion collision frequency. W = Vv + (Vv)^ — 2/31 (V • v) is 
the rate of strain tensor. 

In a dilute astrophysical plasma heat flux Q is dominantly along the magnetic fleld 
hues. The parallel heat flux is given by 



Q = -Xcb(b-V)r, (8) 



where xc is the Coulomb conductivity given by Spitzer (1962) as xc — 6 x 10 '^T^/^ 
ergscm~^K~^ . 
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2.3. LineEirized Expressions for Perturbed Quantities 

We apply a standart Wentzel-Kramers-Brillouin (WKB) perturbation analysis on the 
equilibruim state. To this analysis, all the variables in the MHD equations are denoted by 
sums of an equilibrium value (denoted with a "0" subscript) and a small perturbed quantity 
(denoted with 6) 



p = Po + dp, 
V = Vo + Sv, 
P^Po + SP, 
n = Ho + 5U, 
B = Bo + 5B, 

Q = Qo + SQ. (9) 

Substituting the formulae in equation (9) into equations (2)-(8) and retaining only 
terms up to linear order in perturbations, the linearized perturbation equations are obtained 
as 

V • (5v = 0, (10) 



+ . w = ^vp - V . in - iv (sp + '-^) + ^-^SB, (11) 



dt p^ p \ Air J Airp 



— - = V X (5v X B) + V X (v X 5B) , (12) 



- 10 - 



2)Ot p 2)P 

The perturbed parallel and gyroviscous components of the stress tensor are 



(13) 



5ir = 0.96— 



(I - 3bb)((5b ■ W ■ b) + (I - 3bb)(b • 5W • b) 
+(I-3bb)(b- W-(5b) 



(14) 



5W 



5b X W • (^I + 3bb) + b X 5W ■ ( I + 3bb 
+b X W ■ 3dbb + b X W • 3Mb 
5b X W • + 3bb) + b X 5W • 
+b X W • 35bb + b X W • 3b5b 



(15) 



The perturbed heat flux is given by 



(16) 



The perturbed unit vector of the magnetic field is given by 6h = 6(B/B) 



5B/B -h{dB/B). 



2.4. The Physical Structure of Modes 

Before we obtain the dispersion relation which includes all anisotropic transport effects 
for general axisymmetric disturbances, let us take a glance at the physical meaning of finite 

Larmor radius effect and modes which emerge in the plasma. We consider the local stability 
of a uniformly rotating dilute plasma included only a vertical magnetic field, B — Bz. We 



ignore parallel viscosity, heat flux and radial stratification. We restrict ourself to plane wave 
perturbations that depend only on i.e., of the form exp(iA:2; — iut). Thus, the radial and 
azimuthal components of the linearized motion equation are obtained as 



iujSvR - {2VL - Wgyro)5v^ 



ikB 



SBr = 0, 



(17) 



iujSv^ + {2Q - 4:Vgyro)SvR 



ikB 



SB^ = 0, 



(18) 



4:71 P 



where Vgyro — klP/AudP is the inverse time scale of the gyroviscous stress. The same 
components of the linearized magnetic induction equation are 



The gyroviscous force introduces a term like a Coriolis term in the equation of motion. 
The dispersion relation is 



where v\ — B'^/Anp is the Alfven velocity, vj) — PQ/2uJciP is the drift velocity. 

FLR stress results from changes in particle drift velocities across a gyro-orbit. This 
stress gives rise to distortions of particle orbits and guiding-center drift. Kaufman (1960) 

presented a detailed discussion of this stress. Because ions and electrons have different 
Larmor radii, they move differently due to FLR effects. This different motion gives rise 



— iuSBji — ikBSvR = 0, 



(19) 



— iuoSB^ — ikB6v^ — 0. 



(20) 




2' 



(21) 
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to charge separation. Then it produces a finite parallel electric field. Physically, the FLR 
effects introduce a drift wave that convects the perturbations along the velocity gradient. 

If the angular velocity and magnetic field vectors are oriented in the same sense, the 
term including Vgyro is positive. Prom equations (17)-(20), one can easily see that the 
induced drift motion is opposite to the Coriolis force. Thus, the dynamical epicycle is slowed 
and magnetic tension force is effectively increased. This, in turn, increases the angular 
momentum transfer. The result is an instability. If the angular velocity and magnetic field 
vectors are counter aligned, the signs of these effects should reverse. 

In the limit Q — > 0, the dispersion relation is obtained as 



4 2 
UJ — UJ 



/ Pmc\ ' 



+ k%\ = 0. (22) 



The solutions of equation (22) give two roots: 



± 



1/2 



(23) 



One of the roots describes pure drift mode at large wavenumbers. For small 
wavenumbers (low frequencies) the other root is 



kVf 



) 



(24) 



corresponding to Alfven waves with the gyroviscous force producing a small 
frequency-sphtting of the Alfven wave (see Figure 1). 



One may rewrite the dispersion relation (21) in the presence of the uniform rotation: 
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Fig. 1. — A wavenumber-frequency diagram in the absence of rotation. Forward and back- 
ward waves correspond to the positive and negative values, respectively. 




k 



Fig. 2. — A wavenumber-frequency diagram in the presence of uniform rotation. Forward 
and backward waves correspond to the positive and negative values, respectively. 
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2 , / 4^1) 
± OJUi 1 1^ 



0, 



(25) 



where drift wave frequency is u'j^ = Pmc(k ■ ft)^ /2ep{fl ■ B), Alfven wave frequency is 
uj\= (k- B/47rp), pure inertial wave frequency is cui — 2{k- Q)/k for general magnetic field 
geometries and wavenumbers (Mofatt 1978). The solutions of the Equation (25) are 



1 / ^ujI\ 1 

OU = ±-UJl ( 1 ^ I ± 



1 - 



D 



+ 4a; 



1/2 



(26) 



For (k • $7) = 0, it follows — — uja- If the first term is greater than the second 
one in the brackets in the equation (26), then it is possible to carry out a Taylor series 
expansion of this equation. One then finds a very clear splitting of the fast and slow wave 
frequencies, 



a;+ = ±(jOi I 1 — 



Aid 



D 



( 



\ 



1 + 



-1 4w?, 

a;n 1 ^ 



(27) 



and 



^_ ^ ± (28) 



respectively. These waves result as a combination of the inertial and Alfven waves. 
They are referred to as magneto coriolis (MC) waves. When the magnetic field and rotation 
axis are aligned, the MC wave which imparts a circularly polarized component to the 
velocity perturbation is in the anticlockwise direction. Therefore Coriolis and Lorentz forces 
are in phase. Resulting force causes inertial acceleration. This mode known as fast MC 
wave (i.e., a;+ mode). Indeed, in the equations (27) and (28) the inertial wave is coupled 



- 15 - 



with the drift mode. So, this mode has a larger frequency compared that of the pure fast 
MC wave. For uj_ mode, the MC wave which imparts a circularly polarized component to 
the velocity perturbation is in the clockwise direction. Hence, Coriolis and Lorentz forces 
are out of phase. Resulting force is weakened. This slow MC mode is sometimes referred to 
as magnetostrophic wave or hydromagnetic-inertial wave (see Figure 2). These waves are 
especially important in the dynamo problem (Moffatt 1978, Achenson & Hide 1973). 

When the first term is smaller than the second one in the brackets in equation (26), 
one obtains modified Alfven waves with frequency splitting after Taylor series expansion: 



A wavenumber-frequency diagram of these waves is given in Figure 3. In general the 
angular velocity and magnetic field vectors will not be parallel and the situation will be 
more complex. But the physical picture remains the same. 

2.5. Dispersion Relation with All Anisotropic Transport Effects 

We now consider the axisymmetric behavior of the instability for helical magnetic field 
and more general wavenumbers. 

2.5.1. Equilibrium state 

We work in a cylindrical coordinates system, {R,4>,z). The plasma is assumed to be 
thermally stratified in the presence of a uniform gravitational field in the radial direction, 
g = — g'R. The weak magnetic field is taken to be B (0, BqcosO, BosinO), where 
9 = ta.n~^{Bz/ B^) which is called "pitch angle". In equilibrium the field lines are assumed 




(29) 
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to be isothermal and therefore heat flux is neghgible. We consider differentiaUy rotating 
plasma with a Keplerian velocity profile, i.e. = RQ{R). For this profile and weak 
magnetic field, the stress tensor is negligibly small in the equilibrium for an arbitrary pitch 
angle (see Appendix B). Hence, plasma is in a hydrostatic equilibrium: 

^ = g + i?Q2. (30) 
Po 

2.5.2. General axisymmetric disturbances 

All the perturbed quantities are assumed to have a space-time dependence exp(ik-r-|-a;i), 
where k = kjiR + k^z. WKB assumption requires kR » 1. The time dependence of the 
perturbations are assumed as exp(a;t). This assumption ensures that all coefficients in the 
dispersion relation are real. We work in the Boussinesq limit. In this limit, pressure changes 
are much smaller than temperature and density changes, i.e., ST — —T{5p/p). 

The above equations may be written in explicit component form. The perturbation 
equation of mass continuity is given by 

kuSvn + kzSvz — 0. (31) 

The radial, azimuthal and the axial components of the linearized momentum 

conservation equation are given by the equations (32), (33) and (34), respectively. 



SpdP SP 1 1 

ujSvr - 2fl5v4, ^7^77 + ikn h - — ikR {B^SB^ + BJB^) - - — ik^BJBu 

oR p A-Kp A-Kp 



^par 



^hif^sin 2eSvR + |a sin 29 6v^ + 2^ sin^ eSv, 



Vgyro 



din R oj kz k. 



2sin^dfR^(|f - l)5vR + {2cose^^ + E)5v, - {2^^s\ne + A)5v^ = 0, (32) 
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^2 



dQ 1 
dlnRu 



cosOSvr + cos9SV(j) + sin 95vz 



gyro 



-2D 



dU 1 



an 1 



+ 



0, (33) 



SP 1 

a;5't;2 + ik\ h - — ikzBsSBs + Kar 

p 47rp 



ct m it a; 



dlnRoj 



0. (34) 



Similarly, the radial, azimuthal and axial components of the linearized magnetic 
induction equation are given by the equations (35), (36) and (37) respectively. 



ojSBji — ikzBzSvji = 0, (35) 
uSB^ - ik^BJv^ - 5Br = 0, (36) 
cuSBz — ikzBzSvz = 0. (37) 

Finally the linearized energy equation is, 

t ■ 2n\ . /3dlnPp-5/3 1 . ^ d\nT\ ^ 

— + Vther sm^ e) - 5vR -% + Vther- sm' 9—— = 0, 38 

p ^ ' \5 dR CO dR J 

where Vpar — 0.96klP/2i'p is the inverse time scale of the dissipation due to parallel 
viscosity, Vgyro — klP/AoudP is the inverse time scale of the gyroviscous stress and 
Vther = '^k^xT/^P is the inverse time scale of the dissipation due to thermal conductivity or 
thermal frequency, is epicyclic frequency. Other constants which depend on pitch angle 
^ are ^ = sine(l - 3cos2e), S = sine(l + Qcos^^ - 3sm^9), C = 2sin2e(-l + 381-0^9), 
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D = 3 sin^ 0cose, E^2 cos ^(1 + 3 sin^ 0), F ^2sm0D- sin 20, G = sin ^(1 + 3 cos 20) and 
H — 3sin^cos^ 0. 



2.5.3. Dispersion relation 

The set of equations (31)- (38) are reduced to form M • 5v = 0, where M is a 3 x 3 
matrix. | M |= gives nontrivial solution. Thus, the dispersion relation is obtained as 



Lo^ + a^u^ + a^uj^ + a2up' + aiu + ao = 0, 



(39) 



where 



a4 



(40) 



as 



^2 gyro 



V^^,S^\2sy^ - Af + V,yro^'{G + if) 



(41) 



7 2 1.2 iJl h,2 



+VtherS 



2klv\ + f + |v2^^fi2(25^2 _ ^)2 + v;,,,(]'(G + H) 



+A^f^Vgyro{2sy' - A) 



(42) 
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~ T^VtherS'^Vpar^cD — 



IdPdlnT 



pdR OR 



^VparS'^k'^V^-j^ + VparVgyro^s'^GQ' -j^ + ■i^Vpar^"2cD + -j^Vgyro^H^DQ' 

+{Vgyrof3yQD{2sy^ - A)Q'!^ 



, (43) 



do — VtherS'^ 



y'Kro'^D\^l'f + klv\ + Vgyro^\G/2 + H) 



X 



fe^ 1 ap ainT 



(44) 



where y — kn/k^, ^' — dTt^ /dlnR, s — sin9, c — cosO, k\ — k\-\- k'^cos^O and 
Vgyro = klPi/AVLojcip. = {3/5p){dP/dR){dlnPp-^/^/dR) is Brunt- Vaisala frequency. 



The above dispersion relation is reduced to several previously obtained relations in the 
appropriate limits. Taking Vpar = Vgyro = Vther = = N"^ = , one recovers the result of 
Balbus & Hawley (1991). Taking Vgyro — Vther — — N'^ — 0, one recovers the result of 
Islam & Balbus (2005). Taking A^^ = and Vther — 0, then one recovers the dispersion 
relation given by Paper 1 (2010). If one sets Vpar = Vgyro — — 0, after some algebra, one 
recovers the dispersion relations given by Balbus (2001). 



2.6. Instability Criterion 

The solution of the Equation (39) gives five modes that exist in weakly magnetized, 
dilute plasmas. By analysis of the Routh-Hurwitz criterion, the instability criterion is given 

by 



G 



ktvl + Vgyro^'i^ + H) 



2 2 kl , 
kz'^A + + Vgyro^' 2 



,G klldP dlnT 
¥^dR OR 



' p,2^ " gyro^^ X^" ) ^ 



(45) 
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This complex criterion is simplified, if one considers vertical magnetic field. In this 
case 9 — 90°, D — 0, H = 0,G — —2 and the instabihty criterion is reduced to 



rn^U A V ni 



2 2 kl dn^ _~ dn^ _klldPd\YiT 



< 0. (46) 



In the absence of gyroviscosity and heat conduction, this criterion is reduced to the 
ideal MRI one; in the absence of only gyroviscosity it is reduced to the ideal MTI one. As 
seen from the instability criterion (46), any dynamic instability appears to be dependent 
on the signs of the angular velocity, temperature gradient and gyroviscous force which are 
the free energy sources. Departures from uniform rotation, and isothermality are indeed 
sources of dynamic instability. Gyroviscous force is coupled with angular velocity gradient. 
If one refers to the definition of Vgyro, one clearly sees that gyroviscous force stabilizes or 
destabilizes the modes depending on the sign of f2 • B. 

We consider the case of astrophysical interest, dfl'^/dlnR < 0. If f2 and are aligned 
in the same direction, i.e. 17 • B > 0, the torque term which is the first factor of inequality 
(46) is positive. The instability condition is determined by the sign of the radial force term 
(second factor) in the inequahty (46), that is, if the second factor is negative then instability 
arises: 

2 2_t7 d^'^ _klldPdhiT _kl dQ^ 

This condition differs from the ideal MRI only by additional terms, i.e., gyroviscous 
and temperature gradient on the left hand side of the inequality (47). To ideal MRI, as 
long as dil'^/dlnR < there will be an instability for the small enough k (Balbus & Hawley 
1998). 
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The gyroviscous force acts in the same direction with the magnetic tension force 
because of f2 • B > 0. Therefore, the instabihty is suppressed because the currents 
produced by E x B drift are being out of phase with the current of MRI eigenmode. In the 
astrophysical situation dP/dR < 0, because of the hydrostatic equihbrium. If one assumes 
that the temperature decreases in the direction of gravity, i.e. dlnT/dR > 0, then the 
dilute plasma is completely stable. If one assumes that the temperature increases in the 
direction of gravity, i.e., dlnT/dR < 0, then the dilute plasma may be unstable only if the 
temperature gradient is very steep. 

But in the situation where ft and B^^ are counter aligned, i.e., O-B < 0, the gyroviscous 
force acts in the opposite direction to the magnetic tension force and enhances instability. 
Thus one may find unstable modes even if dlnT/dR > 0, this is because the currents 
produced by E x B drift are being in phase with the current of MRI eigenmode. If one 
assumes that the temperature increases in the direction of gravity, i.e., dlnT/dR < 0, then 
MTI arises. In this case, GvMRI and MTI have equal weights for dominance. It is expected 
that, especially at higher k values, former is the dominant one. In many astrophysical 
plasmas, like cooling white dwarfs and neutron stars, hot accretion flows on compact 
objects release the gravitational potential energy and causes dlnT/dR < 0. But in the 
plasmas having the temperature profile which there are no unstable mode to ideal MTI, i. 
e., dlnT/dR > 0, for example in the cooling flow clusters, GvMRI may operate. 

By referring to the general instability criterion (45), one can argue that perturbations 
with kfi always stabilizes because of the last term in the inequality (45). Interpretation of 
the criterion becomes difficult if we take into account the pitch angles different from 90°. In 
that case, the signs and the ratios of the different terms come into play. Therefore, it is 
more instructive to look at the numerical solutions of the dispersion relation. 
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3. Numerical Solutions of the Dimensionless Dispersion Relation 

The dimensionless dispersion relation is obtained that all the terms of Equation (39) is 
divided by Q^: 

7^ + ^47^ + ^'37^ + + hi + &o = 0, (48) 



where 



h = Q^lr^s'^ + K^^Xs', (49) 



+ |(^;..o)^^^(2V - Ar + + QV-rVr,erX's^§, (50) 



+Y2Kro^yXD 



dlnR 



+ {V-rx'3yD{2sy' - A) 



dlnR F 



2X + + g(V;';.J^X2(2V - Af + V;';.oXf^(G + H) 



+4ftl>-,,(2.y2 _ ^) 



(51) 



X I 



+ 



kldXnO? 
¥ dlnR 



A;2 



X I 



^2 r/1n02 



dlnQ2 /G 



^y'^dlnR \ 2 



+ H 
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61^" 



dliiR fc2 



^ fe2 ^par dlnR 



ther"- 



fc^ gyro 



^ fc2 ^'par ^'sj/ro'J'^ dlni? fc2 ^ k'^ ^ par dlnR 



(52) 



bo = VJlXs 



I * nil 



dln02 



{G/2 + ii-) 



-V- I T/n dlnlJf^G I fcidlnnf_ _ fe| p rp 



^k-^y y^gyro) \ d\nR ) 



,(53) 



where 7 = cj/fl, X = kyjn\ V^^,^ = Vgy^o/X, V^,, = V^ar/^X, V^^, = Vt.^r/^X, 
N = N/n, k = K/n, Pm = (3/5M2)((ilnP/(ilni?) and T = {dhiT/dhiR). M, = v^^/cs is 
the Mach number. 

Dilute plasma condition can be expressed as e = uJci^i ^ 1. From this condition one 
can easily derives the inequality V^"^ » ^gyro- Therefore, we assume V^"^ = 1000. We 
consider the convectively stable plasmas, i.e., > and suppose Pm — —1- Figure 4-8 
show numerical solutions of the dimensionless dispersion relation (Equation 48) under the 
particular assumptions. All the figures are drawn as dimensionless growth rate versus 
wavenumber and 6 for the Keplerian rotational profile. 

Figure 4 shows that there is an instability in small wavenumbers for all pitch angle 
values when only heat conduction is considered. Maximum growth rate of the instability 
is smaller than growth rate of ideal MRI (0.75fi). Indeed, this instability is a MTI which 
emerges in the presence of the helical magnetic geometry. Although the plasma with a 
temperature gradient increasing outward is stable to the ideal MTI, it turns out to be 
unstable when threaded by a helical magnetic field. 

When the angular velocity vector and the component of the magnetic field are 
parallel (f2 tt '^z), gyroviscosity assumes positive values. Figure 5a is drawn only for the 
^gyro — ^ ^^"^ = for the case without heat conduction. This instability is a gyroviscous 
instability which mentioned in Paper I. The mode with small wavenumber < 0.5 is unstable 
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Fig. 3. — A wavenumber-frequency diagram in the presence of uniform rotation. Forward 
and backward waves correspond to the positive and negative values, respectively. 
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Fig. 4. — The growth rates of MTI for helical magnetic geometry. This figure is drawn for 
kR = 0. Instability occurs only at small wavenumbers for all the pitch angles. 
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Fig. 5. — The growth rates of GvMRI. Figure 5a is drawn for = 0. Figure 5b shows growth 
rates versus pitch angle and k^/k^ for X = (15/16)^/^. The mode with any wavenumber is 
unstable for the pitch angles 6 < 50°. 




Fig. 6. — The growth rates of the instability in the presence of anisotropic transport effects. 
Figure 6a is drawn for kji = 0. Mode with large wavenumbers is unstable for the pitch 
angles 10° < 9 < 50°. Figure 6b shows growth rates versus pitch angle and k^/k^ for 
X = (15/16)^/^. The radial wavenumber reduces the growth rates. 
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for all the values of pitch angles. But, for the all wavenumbers, the unstable mode emerges 
only when 9 < 50°. Maximum growth rate of instability is about 3Q. Figure 4 and 5a 
clearly show that GvMRI is more powerful and greater an instability than magnetothermal 
one. 

Figure 5b demonstrates dimensionless growth rates versus pitch angle and kn/kz- In 
this figure, normalized wavenumber of the fastest growing mode of MRI is adopted as 
X — (15/16)^/^. For the pitch angles 9 < 50°, the mode with any wavenumber is unstable, 
but even if maximum growth rate is smaller than the case when kn = 0, it is still greater 
than its correspondent in the ideal MRI case. 

The combined effects of the gyroviscosity and the heat conduction on the instability 
are seen in the Figure 6a. The cases of fl tt and {dlnT / din R) > are considered 
together. While the mode with small wavenumbers (~< 0.5) for all the pitch angle values 
is stable, the one with large wavenumbers becomes unstable only when the pitch angle 
is within the range of 10° < 9 < 50°. Comparison of Figures 5a and 6a shows that the 
temperature gradient acts to reduce the growth rate of the instability. Also, it supresses the 
instability which appears at small wavenumbers (~< 0.5) and small pitch angles {9 < 10°). 

Figure 6b shows the growth rates versus pitch angle and k^/kz for X — (15/16)^/^. 
The radial wavenumber of the unstable mode reduces the growth rates. 

When the angular velocity vector and the component of the magnetic field are 
antiparallel (fi B^), gyroviscosity assumes negative values. For the situation Vgy^^ < 0, 
an interesting situation arises. Instability sets in for the larger pitch angles [9 = 60° — 90°). 
This figure is not included here for spatial economy, because it looks like mirror image of 
the one with Q tt (see Fig. 6a), but the growth rates of instability are larger, i.e. ~ 3fl. 
When we increase the value of V^y^.^, irrespective of its sign, the maximum growth rates of 
the instability become larger, for example, for the VIl^^ = 5, the maximum growth rate is 
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2 90 

Fig. 7. — The growth rates of the instabihty in the presence of the steep temperature 
gradient. Instabihty with higher growth rates arises for 6 < 40°. Besides, there is another 
but rather weak unstable region reveahng itself as a mild crest for 6 > 40°. 




Fig. 8. — The growth rates of the instability for the situation of ujcond ^ ^dyn ^ k ■ v^. 
Figure 8a and 8b are drawn for = 0. The instability is supressed when T = 5 but when 
T = —5 almost at all the wavenumbers and pitch angles instability arises. 
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~ en. 

Now, let us investigate the situation wherein the temperature gradient decreases 
outward. Instabihty behaves in the same manner, that is, the appearance of figures are the 
same but the maximum growth rate is comparatively smaller. 

If temperature gradient assumes a value of 5, and V^^.^ > then instability arises for 
9 < 40°. Besides, there is another but rather weak unstable region revealing itself as a 
mild crest for 9 > 40°. For kvA/n = 2, the growth rate is 0.45 at the pitch angle 50° and 
0.13 at 80°. The growth rates have bigger values at smaller wavenumbers (see Figure 7). 
If temperature gradient assumes a value of —5, instabihty is seen for 9 < 40°, but there 
is no mild crest for 9 > 40°. Growth rates remains as the same. When V^j.^ < and 
dlnT/dlnR — 5, there is instability with the same growth rates at pitch angles 9 > 60°. 
For 9 < 60° again a mild crest appears as a second unstable region. 

In the weak magnetic field limit, cUcond ^ ^dyn ^ k • ordering holds true. Here 
i^dyn ~ {g/Hy^"^ is the local dynamical frequency (Quataert 2008). If t^j^er — Vther/XO, — 
{2klxT/5P)/{kvA) = cOcond/kvA = 1000 and = 1, dlnT/dlnR = 1, then the figure of 
the instability looks exactly like the Figure 5a, that is, only the mode with wavenumbers 
smaller than 0.5 is unstable for all the possible pitch angles, and can not alter the maximum 
growth rate (see Fig. 5a). If this situation is compared with the Figure 6a which is assumed 
^ther = 1 we sec that only the mode with very small wavenumbers is unstable for very large 
characteristic conduction frequency. 

Steep temperature gradient and very large characteristic frequency about conduction 
assumption reveals a very interesting situation. While instability is supressed for T = 5, 
when T — —5, it arises almost at all the wavenumbers and pitch angles (see Figs. 8a and 
8b). 
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3.1. Summary 

We have deduced the effects of various parameters on the behavior of instabihty from 
numerical solutions of the dimensionless dispersion relation. These are 



1. Maximum growth rate of the GvMRI (2.5 — 30) is greater than the growth rates of 
ideal MRI and MTI. 

2. The growth rate of unstable mode depends sensitively on the pitch angle and the 
gyroviscosity parameter Vgyro/^X = V^yro- 

3. Gyroviscosity parameter is positive or negative if angular velocity vector and 
component of the magnetic field are oriented in the same or opposite sense. For these 
situations instability regions in the figures are seen as the mirror images of each other. 
In other words, the instabihty sets in at the smaller pitch angles {9 < 50°) for the 
positive gyroviscosity parameter, but when this parameter is negative instability sets 
in at the larger pitch angles (50° < 9 < 90°). 

4. To ideal MTI, plasma is unstable in which the temperature increases in the direction 
of gravity, i.e. dlnT/dlnR < 0. Whereas, we showed that there is at least one 
unstable mode for a plasma with dlnT/dlnR < or dhiT/dhiR > in the presence 
of helical magnetic field. Instability occurs due to gyroviscous force. 

5. What is the effect of temperature gradient on the GvMRI? First of all, the maximum 
growth rate is reduced. Second, the mode with the very small wavenumbers {k < 0.5) 
remain formally stable (see Fig 5a and 6a). 

6. Only in the presence of the steep temperature gradient, instability occurs at almost 
all the wavenumbers and pitch angles (see Fig. 7). 
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7. In the presence of the steep temperature gradient and very large characteristic 
frequency about conduction, the temperature gradient term stabihzes or destabihzes 
depending on whether it is positive or negative (see Fig. 8). 

4. DISCUSSION AND CONCLUSION 

In a dilute plasma, ion cyclotron frequency greatly exceeds ion-ion collision frequency 
and also electron mean free path is much larger than the gyroradius. Therefore, transport 
of momentum and heat by viscosity and thermal conduction is highly anisotropic with 
respect to the magnetic field orientation. In this regime, for a more accurate plasma 
model the anisotropic transport terms, i.e. parallel viscosity, gyroviscosity and thermal 
conduction, must be taken into consideration in the MHD equations. Many astrophysical 
plasmas display characteristics of the dilute plasma. For example, physical parameters 
and the conditions of the intracluster medium (ICM) of the galaxy clusters are being 
revealed by telescopes with high resolving powers. Chandra X-ray Observatory measured 
X-ray luminosity (10^^ — 10^^ erg/s) emitted by the hot plasma in the ICM and based on 
this measurement the density distribution of ICM as a function of radius is determined 
(Parrish, Stone and Lemaster 2008). Peterson and Fabian (2006) reports that the typical 
densities are in the range of 10~^ — 10"^ cm~^ ; temperatures are 1 — 15 keV. CariUi and 
Taylor (2002) estimated the magnetic field strength in the center of ICM about 1 — 19 /xG 
and 0.1 — 1.0 /iG at the radius of 1 Mpc. With the above quoted values plasma beta is 
P = SnP/B^ ~ 200 - 2000. This value implies that the ICM plasma is dilute and the mean 
free path of electrons is much longer than their gyroradius (Narayan & Medvedev 2001). 

Under these physical conditions, dilute and hot plasma in a differentially rotating disc 
is open to a blend of instabihties hke, MRI, MTI and GvMRI. We proceed the analysis 
from simple to complex. Let us start with a useful mechanical model developed by Balbus 
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& Hawley (1992) for the MRI. This model consists two fluid elements which are tethered 
each other with vertical magnetic field. They are also embedded in a radial angular velocity 
gradient, so that the element orbiting at smaller radius rotates rapidly than the other 
element orbiting at larger radius. We suppose that these elements are residing at different 
vertical locations, but at the same radial locations initially. When these elements are 
radially displaced, the magnetic field will force them to return at their original locations. 
The outward element acquires angular momentum because it has a smaller velocity in its 
new radial location. But the inward element loses angular momentum, because it has a 
greater velocity in its new radial location. When the field fines become more stretched, since 
the inner element continues to lose its angular momentum, it will fall farther inward; the 
other element moves farther out and gains the higher angular momentum. Thus, process 
runs away and instability occurs. This is a classical MRI picture (Balbus & Hawley 1992, 
1998). The ions (because ion viscosity is higher than that of electrons in the dilute plasma) 
in the fiuid elements gyrate around the magnetic field lines and will be under the infiuence 
of a spatially varying electric field which arises from FLR effects. In the shear fiow there is 
no rest frame. Therefore the elements in different locations are exposed to different electric 
fields (E = —V X B/c = —Rfl{R)BQsm9Ii). Along a gyro-orbit, the length scale of electric 
field is comparable to length scale of velocity gradient (WiUiams & Jokipii 1991). If the 
initially positions of elements is selected as a rest frame, the fiuid velocity increases towards 
smaller R values. So, the inward clement will see larger electric field because it has greater 
velocity according to the outward element. Since the drift is controlled by the magnitude 
of the electric field, the fiuid element with a larger relative velocity will drift more rapidly. 
Because the magnetic field which tethers fiuid elements to each other acts as spring-like 
force, the rapidly increasing element separation gives rise to growing spring tension. Thus, 
process runs away and instability occurs quickly. This is a GvMRI picture. 

The modified Hill equations by the inclusion of gyroviscosity, thermal conduction and 
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parallel viscosity can help one to get a better physical understanding of instability. In the 
absence of all three dynamical effects, one recovers the original set of equations describing 
the MRI (Balbus & Hawley 1992, 1998). In equations (54) and (55) below, and are 
the radial and azimuthal displacements of fluid elements. 



where all the abbreviations are the same as above. 

The right hand sides (r.h.s.) of the equations (54) and (55) represent the torque apphed 
to the fluid elements which correspond to "spring" constants in radial and azimuthal 
directions. Comparison of r.h.s. of the equations (54) and (55) show that the gyroviscosity 
couples to the differential rotational, the thermal conduction couples to the Coriolis force, 
the radial gradients of the temperature and the pressure. These complex couplings make 
the roles of the anisotropic forces intangible. We solve equation (55) for and then 
substitute it into equation (54) to find the acceleration of the perturbed fluid element in 
the radial direction. In order the hot, dilute and differentially rotating disc to be unstable 
there should be a net outward acceleration: 



+ {N^cu - n^Vther sin^ OPmT) {cu + Vther sin^ e) 




(54) 




(55) 



{kzVA) in 
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^ + Vther sin 9 



in (56) 



where F — 



Accelaration given by the equation (56) can be written as a sum of three terms 
respectively, i.e., d^^n/dt'^ ^ ar ^ aM + aH,th + 0'M,gv ■ 

Om, given by the equation (56) is the acceleration term arising from the magnetic 
tension force which is always stabilizing. The second acceleration term aH,th is due to the 
radial bouyancy force and anisotropic thermal conduction. And it is a pure hydrodynamic 
term. This term is related to the convective instability in which the source of free energy 
is the temperature gradient. In the absence of the viscous force, when the temperature 
increases in the direction of gravity, this term has a destabilizing effect and when the 
temperature decreases in the direction of gravity, it has a stabilizing effect. The third 
acceleration term aM,gv is due to the radial gyroviscous force. So, it is related to the 
GvMRI. This term depending on the pitch angle, the wavenumbers, gyroviscous force and 
the parallel viscosity may be either positive or negative sign which corresponds to the 
destabilizing and stabilizing effect, respectively. 

The relative importance of the two terms which can be stabilizing or destabilizing 
are depicted in Figure 9. Figure is drawn as the logaritmic normalized wavenumber 
versus logaritmic ratio of accelerations log{aM,gv/ 0'H,th) ^oi convectivcly stable plasma, i.e., 
N'^ — 1. If aM,gv > 0'H,th-i the all unstable modes are GvMRI modes. If aM,gv < 0'H,th, the 
unstable modes are called as thermal modes. As shown in the Figure 9, MTI by driven 
heat conduction is dominant in the very small wavenumbers only for the 6 = 90°, i.e., 
in considering only vertical magnetic field. In all the possible situations throughout the 
whole wavelength range, acceleration term aM,gv due to the gyroviscous force is dominant. 
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Fig. 9. — Figure is drawn as the logaritmic normalized wavenumber versus logaritmic ratio 
of accelerations logiau^gvl (^H,th) for convectively stable plasma, i.e. A^^ = 1, in the presence 
of the three dynamical forces, i.e., parallel viscosity, thermal conduction and gyroviscosity. 
Throughout the whole wavelength range, acceleration term aM,gv which due to gyroviscous 
force is dominant. Also, this acceleration increases rapidly in the large wavenumbers. 



-35 - 



Moreover, this acceleration increases rapidly in the large wavenumbers. {aH,th)CR term may 
be positive for the given values. Still, this term does not make a contribution to instability 
even if V^l^^ = 1000. 

The above analysis shows that, regardless of the sign of the temperature gradient, a 
weakly magnetized and convectively stable dilute plasma harbouring the combined effects of 
gyroviscosity, parallel viscosity and thermal conduction is unstable due to the gyroviscous 
force. Extension of the unstable regions and the growth rates of GvMRI depends sensitively 
on the pitch angle and the gyroviscosity parameter. When the angular velocity vector and 
component of the magnetic field are parallel (O j"!" B^), gyroviscosity assumes positive 
values. When the angular velocity vector and component of the magnetic field are 
antiparallel {fl ti B^), gyroviscosity assumes negative values. Former, while instabihty 
sets in for the smaller pitch angles (6* < 50°), later it sets in for the larger pitch angles 
{9 — 60° — 90°). Maximum growth rates of the unstable mode is approximately 2.5Q and 3Q 
respectively. These values are higher than the ones of magnetorotational instability (0.75Q) 
and the ones of magnetothermal instability (0.5f2, see Figure 4). If the kn wavenumber is 
taken into consideration, the instability region becomes narrower and the growth rate of 
the instability is reduced. 

Briefly, GvMRI is the dominant one in the blend of all three, that is pure MRI, MTI, 
gyroviscous modified MRI. Because the magnitude and the direction of the temperature 
gradient do not exert great effect on the unstable mode, this instability may work at all the 
astrophysical media including dilute plasma. 

I thank E. Rennan Pekiinlii for helpful discussions. I am particularly grateful to 

the referee for constructive comments which led to a substantial improvement of this 
manuscript. This research was supported by the Scientific & Technological Research 
Council of Turkey (TUBiTAK). 
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A. Two-fluid Equations 

The two-fluid equations together with the Faraday and Ampere laws, respectively, are 
jiven by (Braginskii 1965): 



= -n,V ■ v„ (Al) 



^'^'"^ = - VP, - V • n, + q,ns { E + ^ ) + m,n,g, (A2) 



^".^ = -P.(V-v,)-V-Q„ (A3) 



— = -cV X E, (A4) 



J = -^V X B = e {Zui^i - rieVe) , (A5) 

where the subscript "s" stands for electrons and ions, m is the mass, n is the number 
density, v is the velocity of the plasma components, P is the scalar pressure, 11 is the 
stress tensor, q is the particle charge, E and B are the electric and magnetic fields, g is the 
gravitational acceleration, Q is the heat flux, c is the speed of the light, e is the electric 
charge (g = e), Z is the charge state and d/dt — d/dt -\-v - V \s di, Lagrangian derivative. 

Most plasmas of interest are electrically neutral over sufficiently long distance and 
timescales. One assumes that the Debye length is smaller than all 

the relevant spatial scales Xd <^ L. Hence ion and electron densities are essentially equal. 
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i.e. He — Zrii, quasi-neutrality. Then the mass continuity equation is written only for ion 
number density: 



drii 



= -n,:V • V,:. 



(A6) 



Let us multiply the equation (A6) with mj. Since p = mirii, the equation (A6) now 



reads 



dp 
'di 



+ pV • V = 0. 



(A7) 



Let us first write E by using from the electron momentum equation: 



E 



(iVe 

-niene—r- - VPg - V • Xlg + meriej 



dt 



Ve X B 



(A8) 



Let us eliminate E (Equation A8) by substituting from the electron momentum 
equation into ion momentum equation: 



dVj 
dt 



TTiini—r- — — VPj — V • Hi + Zen,, 



1 



dWf. 

-mjle—r- - VPe - V • lie + "le^^eg 

dt 



—Zerii 



X B 



+ Zen, 



Vi X B 



+ niiUig,. (A9) 



Let us put p — niiUi, P — Pe + Pi, U — He + 11% and assume that me/rrii ~ 0, because 
rrie <^ rrii. With the substitutions and using Ampere laws, the ion momentum equation 

becomes 



rfv JxB 

p— = -VP - V • n + + pg. 

dt c 



(AlO) 
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Let us substitute E (Equation A8) into the magnetic induction equation as given (A4): 



dB _ 



crrieZ ( dve\ 1 1 „ „ „ 

g - VPe V • He - Ve X B 

enpUii \ at J eUp eUp 



(All) 



Since it is assumed that Ue — Zrii then J = (c/47r)V x B = eue (v^ — Ve) may be 
written in this form. Prom this equation one obtains Ve = Vj — J/eUe- After substituting Vg 
and using rrie/mi ~ 0, one obtains magnetic induction equation as given below: 

^ = Vx(vixB - —J X B + — VPe + —V ■ n J . (A12) 
at \ eUe erip erip J 

The second term on the right-hand side of the equation (A12) is Hall effect. This term 
is negligible because in the present investigation /3 ^ 1 limit is considered (/3 is the ratio 
of the gas pressure to the magnetic pressure) . The third term on the right-hand side of the 
equation (A12) is thermodiffusion term. The ions carry the most of the momentum due to 
their higher masses. Therefore, the third and the last term on the right-hand side of the 
equation (A12) are negligible. So that resulting magnetic induction equation is given by 

^ = V X (v X B) . (A13) 

If the ions and electrons are in thermal equilibrium, then Tg T. Let us add ion 

energy equation and electron energy equation: 



(JjT (iiT 2 2 J 2 

{Tie + ni)—={Z+ l)n,— = --V • Vi(Pe + P^) + gPeV " — " 3 V " Q- (A14) 

One assumes an ideal gas equation of state. Therefore one may write P — (rig -|- ni)T. 
One obtains the resulting energy equation after some algebraic operations: 
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dP 5^,_ , 2_ ^ 
_ + _P(V.v) = -3V.Q. 



(A15) 



B. The Negligible Effect of Stress Tensor on The Equilibrium State 

On the equilibrium state, the contribution from parallel viscosity components of stress 
tensor is 



(Bl) 



The rate of strain tensor has two components in the equilibrium state, VFr,/, = W^j^R — 
dQ/dlnR. The unit vector along the magnetic field is b = ^b^ + zb^. For this case. 



(b • W • b) = (^Rb^W4,R^ ■ (J>b^ + 56,) = 0. 



(B2) 



Accordingly, there is no contribution of the parallel viscosity to the equilibrium state. 

In the equilibrium state, the contribution from gyroviscosity components of stress 
tensor is 



4:U!ci dR 



b X W • (l + 3bb) + 



b X W 



I + 3bb j 



Since b x W • 3bb = 0, the gyroviscosity component is given by 



(B3) 



v-nr 



P - d 

^' . ( b X W • I + 



Aud dR 



b X W-I 



(B4) 



Using vector and dyadic relationships, one obtains 



V ■ ng'^ 



-sin9 



d 



dR \ d]nR 



+ 



2 dQ 
RdlnR 



The hydrostatic equihbrium equation (30) should be written as 



^^g + Rn' + ^sinO 



d 



dVL 



dR \dlnR 



2 dn 

RdlnR 



Now, Equation (B6) may be rewritten for the gyroviscosity parametre V^'^^o 

QPi/4uJciPVA- 



VP 

Po 



1 d 



dn 



QdR \ dlnR 



+ 



1 2 dfl 
QRdlnR 



The first term in the bracket in the equation (B7) is 



I = 



1 d 



dn 



ndR XdhiR 



2n d 



2n'^dR \d\nR 



dn 



2Q2 



dR \ dlnR 
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dn dn 

'dRdlnR 



because 



-^(2n^ 

dR V dlnR 



2n 



d 



dn 



+ 2- 



dn dn 



dR \d\nRj dRdlnR' 
After some algebraic manipulation Equation (B8) is 



then, 



2Q2 



d / 2(^lnn2 
dR [ dlnR 



-n 



dndinn"^ 

dR dlnR 



2Q2 



dndlnn^ . d 

n h n — 

dR dlnR dR 



dlnn^ 
dlnR 
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then, 



1 dlnn'^dlnn^ d /d\nn^ 
iR dlnR dlnR ^ dR \ dlnR 



(B12) 



The second term in the bracket in the equation (B7)is 



_ 1 2 dn _ n 2 dn _ i dinn^ 

~ QRdlnR ~ Q^RdlnR ~ R dlnR ' ^ ' 

Now, substituting the Equations (B12) and (B13) into the equihbrium state equation, 
we find 



VPo 



g + Rn^ + V^-^^2v\sine 



1 dlnn^dlnn'^ d fdinn^ 

+ 1^ 



m dhiR dlnR OR \ dlnR 



+ 



1 dlnQ^ 
R dlnR 



■ (B14) 



In a Keplerian disc, din fl^ /dlnR = —3. Thus, the equihbrium state equation is 
obtained as 



Y^^S + Rn^-^-V;^J-As^ne. (B15) 

One may consider v\/ R oc B^/R. In a dilute plasma, especially in the intracluster 
medium of galaxy clusters (see Chapter 4), the magnetic field is extremely weak and R is 
relatively very large. Therefore, the contribution of the stress tensor to the equilibrium 
state is negligibly small. And the equilibrium state is given by 



(B16) 
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